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1.  INTRODUCTION 


High-speed  flows  that  have  a  small  lateral  asymmetry  occur  in  rocket  motors 
that  have  canted  nozzles  (to  avoid  impingement  of  the  plume  on  adjacent  struc¬ 
tures)  or  nozzles  that  can  be  vectored.  Interest  in  this  type  of  flow  has  been 
renewed  recently,  since  it  is  thought  to  play  a  role  in  an  instability  observed 
in  certain  spin-stabilized,  solid-propellant,  upper-stage  rockets.  In  this 
case,  the  flow  at  the  nozzle  entrance  is  oblique  to  the  nozzle  axis  as  a  conse¬ 
quence  of  the  Coriolis  force  acting  on  the  combustion  gas  in  a  precessing  and 
nutating  rocket  (Ref.  1,2). 

Nonaxi symmetric  flows  through  nozzles  have  also  been  considered  by  Walters, 
(Ref.  3)  both  experimentally  and  theoretically.  However,  in  that  case,  the 
nozzles  themselves  were  nonaxisymmetric ,  having  an  obliquely  machined  throat 
section.  Other  authors  have  studied  nonaxi symne trie  flows  in  nozzles,  but 
generally  they  have  confined  themselves  to  the  supersonic  part  of  the  flow, 
computed  by  the  method  of  characteristics,  starting  from  an  assumed  inclined 
sonic  surface  at  the  nozzle  throat  (Refs.  4,5).  None  of  these  studies  is 
directly  applicable  to  the  case  encountered  in  canted  nozzles  or  in  unstable 
spinning  and  precessing  rocket  vehicles. 

In  these  applications,  the  flow — as  a  correction  to  the  axisymmetric  flow — 
is  sufficiently  closely  approximated  by  assuming  a  steady-state,  inviscous  and 
adiabatic  flow  of  a  perfect  gas.  Based  on  the  assumption  of  constant  reservoir 
conditions  upstream  of  the  nozzle,  the  flow  is  isentropic  and  isoenergetic 
(constant  total  enthalpy)  everywhere.  Analogously  to  the  well-known,  elemen¬ 
tary  theory  of  quasi-one-dimensional  flow,  the  nozzle  cross  section  is  assumed 
to  vary  only  slowly  with  axial  distance.  Since  the  lateral  asymmetry  caused  by 
the  angle  of  attack  or  transverse  pressure  gradient  of  the  entrance  flow  is 
also  assumed  small,  the  effects  due  to  the  changing  cross  section  and  those 
due  to  the  lateral  asymmetry  can  both  be  treated  as  small  perturbations 
superposed  on  the  one-dimensional  zero-order  flow. 
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and  similarly  for  the  density  p,  enthalpy  h,  and  the  Cartesian  velocity 
components  u,v,w;  hence,  for  instance,  for  u  we  write 

u  =  uq(z)  +  eu^  +  s.f.p.  (lb) 

The  symmetric  flow  perturbation  (s.f.p.)  terras  are  written  informally  merely  as 
a  reminder;  it  will  be  clear  from  the  subsequent  development  that,  by  reason 
of  symmetry,  they  do  not  contribute  to  any  of  the  integrals  defined  below. 

The  first-order  asymmetric  terms  p^,h^,w^  can  all  be  expressed  in  terms  of 

p  .  Thus,  from  the  assumption  of  a  perfect  gas  and  constant  entropy,  p/p  (z)  = 
1  /  Y  ^ 

tp/pQ(z)]  1  where  y  is  the  ratio  of  the  specific  heats.  Since  this  relation 
applies  separately  to  the  symmetric  and  the  asymmetric  case  and  superposition 
holds,  it  follows,  after  neglecting  second-  and  higher-order  terms,  that 

pl/po(z)  =  y  pl/po(z)  (2a) 

and  similarly 

h1/hQ(z)  =  ^  ~  1  p1/pQ(z)  (2b) 

From  conservation  of  energy,  h  +  l/2(u^  +  v^  +  w^)  =  const.,  and  neglecting 

2 

again  second-order  terms,  it  follows  that  w^/wq  =  -h]/w0  •  Introducing  the 
Mach  number  M  (z)  of  the  quasi-one-dimens ional  flow,  M  =  w  /  (y  -  l)h  , 
we  can  express  the  asymmetric  perturbation  of  the  axial  velocity  component  by 

w1/wq(z)  *  -  [yM2(z)]_1  px/po(z)  (2c) 
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2.  CONSERVATION  EQUATIONS 


We  consider  a  control  volume  bounded  by  transverse  planes  at  z  and  z  +  dz 
and  by  the  nozzle  walls.  We  designate  by  eP^x(z)  the  transverse  momentum 
(which  is  in  the  x-direction)  carried  by  the  flow  through  the  transverse 
plane,  per  unit  time.  In  the  case  of  the  rectangular  nozzle,  is  taken 

per  unit  width  in  the  y— direction.  From  conservation  of  momentum  to  the 
lowest  significant  order 


dP 


lx 


dz 


<Pl)xrfl(z)  *  (pl*x— R(z) 


if  a  =  1 


(3a) 


dP 


2ir 


lx 


dz 


=  -R(z)  J  (p1)rssR(z)  cos*  d* 
*=0 


if  cr  =  2 


(3b) 


for  the  rectangular  (o  *  1)  and  axisymmetric  (o  =  2)  nozzle,  respectively. 

The  transverse  angular  momentum  about  the  origin  of  the  coordinate  system 
is  in  the  y— direction.  The  amount  carried  through  a  plane  normal  to  the 
nozzle  axis,  per  unit  time  (and  per  unit  width  in  the  case  of  the  rectangular 
nozzle),  is  designated  by  eL.  (z).  When  second—  and  higher— order  terms 

L*  2 

are  dropped,  the  momentum  flux  pw  in  the  z-direction  becomes 


2  2  2 
Pw  =  P  **  +  2ep  w  w,  +  cw  p,  +  s.f.p. 

oo  ool  orl  r 


By  symmetry,  the  zero  and  symmetric  first-order  terms  in  pw2  make  no 
contribution  to  the  integration  over  x,  so  that  for  the  rectangular  nozzle 

R(z) 

Lly  =  zPlx  ~f  (2poVl  +  Vl)xdx  if  o  =  1  (4a) 

Xa-R(z) 


The  first  term  on  the  right  comes  from  the  u  component  of  the  velocity, 
integrated  over  the  cross  section;  the  second  term  is  from  w. 
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By  essentially  the  same  argument,  for  the  axisymmetric  nozzle 


R(z)  2ir 

L.  =  zP,  -  /  /  (2p  w  w.  +  w2p. )r2cos$  dr  d$ 

*!«  001  Oi 


ly  lx 


r=0  $=0 


Conservation  of  angular  momentum  then  requires  that 


dL. 

dz 


R(z) 


_d 

dz 


/ 


dP 


pjXdx  +  ( z  +  R- 


lx 


dz 


x=-R(z ) 


if  a  =  2  (4b) 


if  o  =  1  (5a) 


R(z) 

if  f 


2ir 


p^r  cos$  dr  d$  + 


z  + 


R 

dz/  dz 


r=0  $=0 


if  a  =  2  (5b) 


where  Eqs.  (3a)  and  (3b),  respectively,  have  been  used  and  where  again  the 
zero-order  and  first-order  symmetric  terms  make  no  contribution  to  the 
integrals.  The  first  term  on  the  right  represents  the  moment  from  forces 
acting  on  the  transverse  planes  bounding  the  control  volume;  the  last  term 
results  from  the  wall  pressure. 


The  term  is  eliminated  from  these  equations  by  differentiating 
Eqs.  (4)  and  substituting  the  result  into  Eqs.  (5).  When  Eqs.  (2)  also  are 
used  to  eliminate  p^  and  w^,  it  follows  that 


R(z) 


x=*-R(z) 


i  v  .  A  D  dR  dPlx 

-  l)p.xdx  +  R  - - - — 

rl  dz  dz 


if  9  *  1 


_d 

dz 


R(z)  2ir 


(M2  -  l)Plr2cos^  dr  d*  +  R^f  -  Plx 


0 


if  a  =  2 
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where  use  has  been  made  of  the  expression  for  the  dynamic  pressure  of  the 
quasi-one-dimensional  flow  (Ref.  6)  (po/2)wq  =  (y/2)p  Mq.o  Differentiating 
again  with  respect  to  z  and  using  Eqs .  (3)  gives  the  following  equations  for 
p^,  where  now  *11  other  first-order  perturbation  terms  have  been  eliminated: 


-.3 - 1 


R(z) 

/ 

p^xdx 

x=-R(z) 

(■*)] 

(P1>: 

R(z) 

r 

2w 

f 

J  . 

r=0 

1  1 

<t>=o 

(■»)] 

2ir 

*/  1 
4>=o 

„dR  _d  /  \  /v 

'  Kdz  dz  'P1  x=R(z)  ~  Vpl;x=-R(z) 


l'x=R(z) 


>l)xs-R(z)  =  0 


if  o  =  1 


2ir 

_dR  d  _  f  .  . 

-  "Si  Tz  *J  (pi>  „ 

L  ♦=<>  r=t 


cos$  d<p 


if  o  =  2 
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3 .  INTEGRAL  METHOD 


Integral  methods  are  typically  based  on  prescribing  for  the  dependent 
variable  a  simple  functional  form  that  may  depend  on  one  or  several  parameters 
and  satisfies  the  boundary  conditions.  The  parameters  then  are  determined 
such  that  the  integral  relations — in  this  case,  the  conservation  equations  for 
the  transverse  momentum  and  angular  momentum  integrated  over  the  nozzle  cross- 
section — are  satisfied.  With  a  judicious  choice  of  the  functional  dependence, 
useful  results,  although  of  limited  accuracy,  can  often  be  obtained. 

In  the  present  case,  we  prescribe  for  the  asymmetric  perturbation  term  p^ 
a  linear  dependence  on  the  transverse  coordinate,  in  the  form  therefore 

P,  gj(z)  x 

iTuT  =  R(z)  (7) 

where  the  nondimensional  coefficient  g^(z)  is  to  be  determined. 

We  designate  by  R*  the  nozzle  half  width  (o  =  1)  or  nozzle  radius 
(a  *  2)  at  the  throat  (a  fictitious  throat  if  the  quasi-one-dimensional  flow 
is  subsonic  throughout)  where  M  =  1,  and  pQ  -  pQ  and  define  the 
nondimensional  quantities 

R*  Po(z) 

S  =  R/R*  ,  C  =  z/R*  ,  kx(0  =  gx(z)  (8) 


Carrying  out  the  integrations  in  Eqs.  (6)  yields  a  combined  relation  for  the 
rectangular  and  axisymmetric  nozzles,  in  the  form  of  a  second-order 
differential  equation  for  k^(£) 


1 

o  +  2 


_d  d£  fO+1  . 

dC  Id;  *  *1 


+  kL  =  0 


(9) 
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Equation  (9)  represents  the  principal  result  of  Sections  1  through  3  and  is 
suitable  for  computer  programming  of  the  direct  problem  (prescribed  zero  order 
Mach  number  as  a  function  of  axial  distance)  or  of  the  inverse  problem 
(prescribed  nozzle  contour).  In  either  case,  it  is  advantageous  to  express 
the  function  £(£)  in  Eq.  (9)  by  the  Mach  number  Mq(£).  From  a  well-known 
formula  for  quasi-one-dimensional  flows  (Ref.  6) 


K  = 


1  2 
M2  Y  +  1 


(Y+1)/(Y-1) 


l/(2c) 


(10) 


As  one  would  expect,  Eq.  (9)  exhibits  the  transition  at  M  *  1  between 
an  essentially  exponential  behavior  of  k^(C)  and,  on  the  other  hand,  a 
wavelike  character  as  the  sign  of  the  coefficient  of  the  highest  (second) 
derivative  changes. 

Boundary  Conditions 

As  is  well  known,  classical  one-dimensional  nozzle  theory  requires  the 
Mach  number  to  be  unity  at  a  throat  where  a  transition  from  subsonic  to 
supersonic  flow  occurs.  An  analogous  condition  occurs  in  the  present  case  and 
is  caused  by  the  vanishing  of  the  coefficient  in  Eq.  (9)  of  the  second  deriv¬ 
ative  of  k^(£)  at  M  ■  1.  Carrying  out  the  differentiations  indicated 
in  Eq.  (9)  and  letting  at  the  throat  M«l,  £*1,  d£/d£  =  0  results 
in  the  condition  for  the  solution  to  be  regular  at  (  *  0 


2 

o  +  2 

This  condition, 
sonic  throat. 


dM2  dk. 

— 2  — L  + 

dC  d(,  + 
therefore. 


relates  the  first  derivative  of  k^ 


at  C  ■  0 


to  k^  at  a 


(lla) 
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A  second  boundary  condition  results  from  prescribing  the  transverse 
pressure  gradient  edp^/dx  at  an  initial  (upstream)  location  £  =  hence  the 

9Pt 


condition 


k,  =  — r 


R* 

* 

P„ 


Ji 

dx 


at  C  =  Ci  (lib) 


In  applications  to  rocket  motors,  the  transverse  pressure  gradient  at  the 
entrance  to  the  nozzle  needs  to  be  determined  from  pressure  matching  with  the 
flow  field  in  the  motor  chamber  adjacent  to  the  nozzle.  Since  the  combustion 
gas  velocity  in  the  chamber  is  typically  far  below  the  speed  of  sound,  the 
calculation  is  simplified  by  the  assumption  of  incompressibility,  although  it 
is  sometimes  complicated  by  the  geometry  of  the  boundaries.  This  latter 
calculation  is  outside  the  scope  of  the  present  report. 


After  k.(£)  has  been  obtained  from  a  solution  of  Eq.  (9),  the  trans¬ 
verse  pressure  gradient  is  obtained  from  Eqs.  (7)  and  (8).  In  turn,  one 
obtains  the  asymmetric  force  eF1  (per  unit  axial  length)  exerted  by  the 
gas  on  the  nozzle  wall.  This  force  is  often  the  quantity  of  greatest  interest 
in  applications.  For  instance,  for  the  axisymmetric  nozzle,  after  integrating 
over  the  polar  angle  $ 

Fx  =  ir(p*/R*)  R2(C)kl(C)  (12) 

Numerical  results  for  the  homogeneous  case,  applied  to  a  parameterized 
family  of  nozzle  contours,  are  discussed  in  Section  5. 
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4.  A  TEST  CASE 


Supersonic  flow,  incident  at  a  small  angle  on  a  rectangular  duct  of 
uniform  width  (Fig.  1),  represents  an  example  in  which  an  exact — at  least  in 
the  sense  of  a  first-order  perturbation  result — solution  is  easily  calculated. 
It  is  of  interest,  therefore,  to  compare,  as  a  test  case,  results  obtained 
from  Eq.  (9)  with  the  exact,  but  more  restricted,  solution  in  this  special 
case.  The  exact  first-order  solution,  which  is  well  known,  is  described  by  a 
periodic  pattern  of  triangular  and  rhombic  regions  bounded  by  Mach  lines,  in 
which  the  flow  properties  are  constant,  but  changing  discontinuously  across 
the  Mach  lines. 


Note:  R  =  half  width;  =  incident  flow  velocity;  X  »  period 


Figure  1.  Supersonic  Flow  with  Small  Angle  of  Incidence  Through 
Rectangular  Duct 
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The  transverse  pressure  gradient,  and  similarly  the  other  perturbation  quanti¬ 
ties,  therefore,  have  a  sinusoidal  dependence  on  the  axial  coordinate,  with 
period 


(integral  method) 


(13b) 


The  integral  method,  therefore,  gives  the  correct  Mach  number  dependence 
although  with  a  multiplier  which  differs  from  the  correct  one  by  approximately 
10X. 
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5.  NOZZLES  WITH  POLYNOMIAL  MACH  NUMBER  DEPENDENCE 


As  an  example  of  the  application  of  Eq.  (9),  we  consider  an  axisymmetric 
Laval  nozzle  for  which  the  square  of  the  unperturbed  (quasi-one-dimensional) 
Mach  number  varies  with  axial  distance  as  a  second-degree  polynomial ,  for 
which  we  choose 


M2  -  1  «  aC(l  -  bC),  a  >  0,  b  >  0  (14) 

o 

where  a  and  b  are  constants.  The  nozzle  contour  is  easily  calculated  from 
Eq.  (10).  For  suitably  chosen  constants,  relation  (14)  results  in  quite 
realistic  nozzle  contours  (Fig.  2). 


r 


Figure  2.  Nozzle  Contours  Satisfying  Eq.  (14),  for  y  *  1*20 
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At  M  =0  the  slope  of  the  nozzle  contour  becomes  infinite.  This  occurs 
o  __  r 

at  C  *  -(2b)-^(/l  +  4b/a  -  1).  Clearly,  the  assumption  of  quasi-one-dimensional 

flow  breaks  down  for  £  approaching  this  value.  At  £  =  (2b)  *,  the  Mach 

number  M  and,  hence,  the  nozzle  radius  each  reach  a  maximum,  which  occurs 
o 

in  the  supersonic  part  of  the  flow. 

The  transverse  pressure  gradient  in  the  form  of  the  function  k^(C)  and 
the  lateral  force  per  unit  axial  length  exerted  on  the  nozzle,  in  the  form  of 
the  nondimens ional  ratio  F^O/Cp***)*  are  computed  for  several  values  of 
the  nozzle  parameters  a  and  b.  They  are  graphed  in  Figures  3  and  4.  (For  the 
range  of  b  that  is  of  practical  interest  for  realistic  nozzle  contours,  the 
curves  in  Figure  3  very  nearly  coincide  and  depend  only  on  the  parameter  a.) 

The  solutions  are  computed  for  a  ratio  of  the  specific  heats  y  =  1.20,  a 
value  that  is  representative  of  many  rocket  motor  combustion  gases,  and  are 
normalized  so  that  k^  ■  1  at  the  nozzle  throat.  It  is  evident  that  the 
largest  transverse  pressure  gradients  and  nozzle  side  forces  occur  in  the 
subsonic  and  transonic  sections  of  the  nozzle  where  most  of  the  readjustment 
of  the  flow  direction  takes  place.  It  is  also  evident,  particularly  for  the 
smaller  values  of  the  parameter  a,  that  the  transverse  pressure  gradient  down¬ 
stream  of  the  throat  at  first  reverses  sign,  an  indication  of  the  reflection 
on  the  nozzle  walls  of  the  Mach  cones  associated  with  the  turning  of  the  flow. 

Table  1  lists  numerical  values  for  the  side  force  on  the  nozzle.  The 
effect  of  Mach  cone  reflections  reversing  the  sign  of  the  force  in  the 
diverging  (supersonic)  part  of  the  nozzle  is  again  evident. 

In  some  applications,  such  as  in  spin-stabilized  solid-propellant  rockets, 
the  condition  at  the  nozzle  exit  plane  is  of  particular  interest,  since  it 
determines  the  angular  momentum  flux  of  the  combustion  gas  and  hence  the 
resulting  moment  about  the  vehicle's  center  of  mass.  Table  2,  therefore, 
lists  the  transverse  force  F1  at  the  exit  plane,  per  unit  axial  length, 
relative  to  the  force  F^  at  the  nozzle  throat,  for  nozzles  with  the 
polynomial  Mach  number  dependence  (14)  and  zero  nozzle  divergence  at  the  exit 
plane. 
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y.d/id 


Figure  3.  Function  k^  (O  Computed  from  Eq.  (9)  and  Normalized 
to  ki  *  1  at  Nozzle  Throat 


Figure  U.  Nondimensional  Ratio  Fi(0/(p*R*)  of  the  Transverse 
Force  F^  on  the  Nozzle  Wall,  per  Unit  Axial  Length 
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Table  1.  Function  Fi(C)/(p*R*)  for  the  Nozzle  Contours  Shown  in 
Figure  2  for  y  =  1.20 


-0.6  +  0.918 

-0.6  +  2.776 

+  3.315 
+  3.161 
+  2.627 
+  2.017 
+  1.636 
0.8  +  0.960 

+  0.562 
+  0.238 
+  0.015 

-  0.163 

1.8  -  0.250 

-  0.319 

-  0.381 

-  0.366 

-  0.318 

-  0.267 

.5  -  0.220 

5.0  -  0.181 

-  0.122 

-  0.086 

8.0  -  0.060 

9.0  -  0.066 

10.0  -  0.036 


1 


a  -  2.00 
b  -  0.020 

a  -  2. SO 
b  -  0.020 

mm 

38 

a  -  2.00 
b  -  0.030 

a  ■  2.50 
b  -  0.030 

a  -  1.50 
b  -  0.050 

a  -  2.50 
b  -  0.050 

a  «  3.50 
b  -  0.050 

4 

0.861 

4- 

0.750 

4 

1.327 

4 

2.717 

♦ 

1.300 

+ 

2.663 

4 

2.810 

♦ 

2.368 

4 

3.305 

4 

2.801 

4 

2.359 

4- 

3.287 

4- 

2.362 

4- 

1.666 

4 

3.161 

4 

3.161 

4 

3.161 

♦ 

3.161 

4 

3.161 

4- 

3.161 

4- 

3.161 

4- 

3.161 

4 

2.879 

4 

3.067 

4 

2.632 

♦ 

2.886 

4 

3.052 

4 

2.660 

4- 

3.062 

4 

3.263 

4 

2.391 

4 

2.616 

+ 

2.022 

4 

2.398 

♦ 

2.625 

4 

2.033 

4 

2.662 

4 

2.869 

4 

1.870 

♦ 

2.118 

4 

1.661 

4 

1.877 

♦ 

2.129 

4 

1.669 

4- 

2.150 

4 

2.365 

4 

1.399 

4 

1.662 

4 

0.962 

4 

1.606 

4- 

1.673 

4- 

0.966 

4* 

1.695 

4 

1.898 

♦ 

1.008 

4 

1.282 

4 

0.561 

4 

1.013 

4 

1.292 

4- 

0.538 

4 

1.316 

4- 

1.538 

4 

0.697 

♦ 

0.980 

4 

0.236 

+  0.699 

♦ 

0.989 

4- 

0.223 

4 

1.007 

4- 

1.260 

4 

0.657 

4 

0.766 

4 

0.007 

+  0.656 

4 

0.753 

- 

0.011 

4- 

0.766 

4- 

1.068 

4 

0.276 

4 

0.568 

- 

0.155 

♦ 

0.271 

4 

0.572 

- 

0.179 

4- 

0.579 

4- 

0.886 

♦ 

0.161 

4 

0.632 

- 

0.265 

♦  0.133 

4 

0.633 

- 

0.295 

4- 

0.635 

4* 

0.762 

4 

0.061 

4 

0.329 

- 

0.335 

4 

0.031 

4 

0.328 

- 

0.370 

4- 

0.325 

4- 

0.667 

- 

0.101 

4 

0.168 

- 

0.601 

- 

0.118 

4 

0.161 

- 

0.662 

4* 

0.165 

4 

0.508 

- 

0.159 

♦ 

0.085 

- 

0.382 

- 

0.180 

4 

0.075 

- 

0.622 

4* 

0.050 

4- 

0.617 

- 

0.178 

4 

0.062 

- 

0.336 

- 

0.200 

♦ 

0.029 

- 

0.366 

- 

0.003 

4- 

0.362 

- 

0.178 

♦ 

0.019 

- 

0.279 

- 

0.202 

4* 

0.005 

- 

0.301 

- 

0.032 

4- 

0.327 

- 

0.171 

+ 

0.006 

- 

0.228 

- 

0.195 

- 

0.009 

- 

0.238 

- 

0.050 

4- 

0.306 

- 

0.162 

- 

0.001 

- 

0.183 

- 

0.185 

- 

0.017 

- 

0.181 

- 

0.061 

4- 

0.289 

- 

0.163 

- 

0.007 

- 

0.116 

- 

0.165 

- 

0.026 

- 

0.088 

- 

0.072 

4* 

0.275 

- 

0.127 

- 

0.009 

- 

0.072 

- 

0.169 

- 

0.026 

- 

0.020 

- 

0.079 

4- 

0  275 

- 

0.116 

- 

0.010 

- 

0.063 

- 

0.137 

- 

0.027 

4- 

0.033 

- 

0.085 

4 

0.287 

- 

0.108 

- 

0.010 

- 

0.023 

- 

0.128 

- 

0.027 

4- 

0.070 

- 

0.093 

4 

0.311 

0.101 

0.010 

0.009 

“ 

0.123 

-* 

0.027 

4 

0.103 

- 

0.103 

4 

0.351 

Table  2. 


Ratio  F^g/F*  of  the  Transverse  Force  F^g  on  the  Nozzle 
Wall,  per  Unit  Axial  Length,  at  the  Nozzle  Exit  Plane, 
Relative  to  the  Corresponding  Force  F*  at  the  Throat 


Fle/F* 

a  =  1.50,  b  -  0.020 

-  2.468  x  10“3 

a  =  2.00,  b  =  0.020 

-  3.811  x  10~2 

a  *  2.50,  b  =  0.020 

-  4.730  x  10-3 

a  =  1.50,  b  =  0.030 

+  9.298  x  10“3 

a  =  2.00,  b  =  0.030 

-  4.230  x  lO-2 

a  »  2.50,  b  =  0.030 

-  1.205  x  lO"2 

a  =*  1.50,  b  *  0.050 

+  3.267  x  10-2 

a  a  2.50,  b  a  0.050 

-  3.279  x  lO"2 

a  =  3.50,  b  =  0.050 

+  1.116  x  10"1 

Note:  For  y  *  1.20,  and  zero  nozzle  divergence 
angle  at  the  exit  plane. 
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a,  b 
F^z) 

81.  kl 
h 

Lly(Z ) 

Mo(z) 

P 

Plx(z) 

R(z) 


u,v,w 


x,y,z;  r,4>,z 
Y 

p 

o 

(  >o 
(  >1 
(  )* 

(  >i 
(  )e 


NOMENCLATURE 
parameters  [Eq.  14)] 

transverse  force  on  nozzle,  per  unit  social  length 

functions  characterizing  transverse  pressure  gradient 
[Eqs.  (7)  and  (8)] 

enthalpy 

transverse  component  of  angular  momentum 
Mach  number  of  zero-order  (unperturbed)  flow 
pressure 

transverse  component  of  momentum 

half-width  of  rectangular  nozzle,  and  radius  of 
axisymmetric  nozzle 

Cartesian  velocity  components 

Cartesian  and  cylindrical  coordinates 

ratio  of  specific  heats 

nondimensional  lengths  [Eq.  (8)] 

density 

o  =  1  for  rectangular  nozzle;  o  =  2  for  axisymmetric  nozzle 

zero-order  term 

first-order  asymmetric  term 

condition  at  Mq  ■  1 

initial  condition 

nozzle  exit  plane  condition 
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